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1. Motivation 

A number of mesons observed in experiments like LHCb or Belle are not well understood. 
Those mesons have masses and quantum numbers, which are not typical for standard quark- 
antiquark states, but indicate an exotic four-quark structure. Prominent examples are the charged 
charmonium-like and bottomonium-like states Z, ± and Z/, + (cf. e.g. [1]). Their masses and decay 
products suggest the presence of a cc or bb pair, respectively. On the other hand their electric charge 
indicates additionally a light quark-antiquark pair ud or du. Those four-quark systems, in the fol¬ 
lowing also referred to as tetraquarks, are expected to be studied in more detail in the near future 
by experimental collaborations. Therefore, a sound theoretical understanding of those systems is 
crucial and of great interest. 

Here we summarize the main results of our recently published work [2], where we have studied 
four-quark systems with two heavy antiquarks bb and two lighter quarks qq using lattice QCD and 
the Born-Oppenheimer approximation. First bb potentials in the presence of lighter quarks qq are 
computed. Then the Schrodinger equation is solved using these potentials, where possibly existing 
bound states indicate stable tetraquarks. Other papers studying the same systems with similar 
methods include [3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14]. 

2. Qualitative discussion of qqbb systems 

At small bb separations the bb interaction is dominated by 1-gluon exchange. For a bound 
state the bb pair must, therefore, be in an attractive color triplet. Due to the Pauli principle and 
because we assume a spatially symmetric x-wave, bb has to form an antisymmetric color-spin-flavor 
combination and, hence, a symmetric spin combination, i.e. bb spin jb = 1. Since the complete 
four-quark system is color neutral, the light quarks qq must be in an antisymmetric color antitriplet. 
Again due to the Pauli principle qq has to form an antisymmetric color-spin-flavor combination 
and, hence, a symmetric spin-flavor combination. Candidates for tetraquarks are, therefore, the 
(spin) scalar isosinglet (i.e. a qq spin singlet j = 0 with antisymmetric flavor, e.g. qq G {{ud — 
du)/s/l , (s^V^ — s^s^)/s/l, (c^M 2 ) — c^c^) / s/l} 1 ) and the (spin) vector isotriplet (i.e. a 
qq spin triplet j = 1 with symmetric flavor, e.g. qq G {uu , {ud + du) / s/2. dd. ss , cc}). The overall 
quantum numbers of a bound qqbb system are I{J P ) = 0(1 + ) for the scalar isosinglet channel and 
I{J P ) G (1(0 + ), 1(1 + ), 1(2 + )} for the vector isotriplet channel. 

At large bb separations the bb interaction is screened by the light quarks qq, i.e. the four 
quarks form a system of two heavy-light mesons. One expects stronger screening for increasing 
quark mass m q , because the wave functions of the corresponding mesons qb are then more compact. 

3. Lattice QCD computation of static antiquark-antiquark potentials 

We extract potentials of two static antiquarks <22 (approximating the two b quarks of the qqbb 
system) in the presence of two light quarks qq from correlation functions 

C{t,r) = {Q\@\t)@{0)\Q.) c< exp(-V(r)f). (3.1) 

t —>°° 

1 To be able to study flavor antisymmetric qq combinations with q = s, we consider two hypothetical degenerate 

flavors with the mass of the 5 quark, T 1 ) and T 2 } and similarly for q = c, F 1 * and A 2 k 
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6 denotes a four-quark creation operator, 

0 = (^T) AB ( < rf) C D(<2c(ri)<?i 1 ) ( r i)) (<2D(r 2 )4 2 , ( r 2)) > r = |i*i — «* 2 |, (3.2) 

where T is an appropriate combination of y matrices accounting for defined quantum numbers light 
quark spin \ j z \, parity P and P x (cf. [8] for details), f G {(1 — 70)75, (1 — 70 ) 7 /} does not affect 
the resulting potential V ( r ), since the static quark spin is irrelevant. 2 ? = 707? denotes the charge 
conjugation matrix. Note that operators like (3.2) generate overlap not only to mesonic molecule 
structures, but also to diquark-antidiquark structures [15, 16]. 

The asymptotic value of a potential and whether it is attractive or repulsive depends on the 
quantum numbers (\j z \.P.P x ) and, hence, on T. In the following we are exclusively interested in at¬ 
tractive potentials between two ground state static-light mesons: the scalar isosinglet corresponding 
to T = (1 + 70)75 and the vector isotriplet corresponding to T = (I + 70)7/. 

Computations have been performed using two ensembles of gauge link configurations gener¬ 
ated by the European Twisted Mass Collaboration (ETMC) with dynamical 11/d quarks. Informa¬ 
tion on these ensembles can be found in Table 1 and [17, 18]. 

/j lattice size /i/ a in fm m n in MeV # configurations 
3.90 24 3 x 48 0.00400 0.079 340 480 

4.35 32 3 x 64 0.00175 0.042 352 100 

Table 1: Ensembles of gauge link configurations (/3: inverse gauge coupling; /.//: bare u/d quark 
mass in lattice units; a: lattice spacing; m n : pion mass). 


4. qqbb tetraquarks in the Born-Oppenheimer approximation 

To determine an analytical expression for the QQ potential or equivalently bb potential, we fit 
the ansatz 



with respect to a, d and Vo to the lattice QCD results obtained in the previous section. The constant 
Vo accounts for twice the mass of the ground state static-light meson. 

We insert the the analytical expression (4.1) in the Schrodinger equation for the radial coordi¬ 
nate of the two b quarks (which we assume to be in an 5-wave), 

{~k^ +m ) R(r) = EsR{r) (42) 

with U(r) = V(r)|y 0=0 and /./ = nip/ 2 and determine the lowest eigenvalue Eg. If Eg < 0, the 
four quarks qqbb can form a tetraquark. If Eg > 0, there is no binding, i.e. the four-quark system 
will always be a system of two unbound B mesons. Notice that this so-called Born-Oppenheimer 
approximation is valid for m q <C «?/,, which is certainly the case for q £ {u. d . .v} and at least 
crudely fulfilled for q = c. 

To quantify the systematic errors of different channels (scalar isosinglet and vector isotriplet, 
different light flavors q G {u. d. s , c}), we perform a large number of fits varying the range of 
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temporal separations t mm < t < t max of the correlation function C(t,r) (cf. eq. (3.1)), at which the 
lattice potential is read off, as well as the range of spatial bb separations r mm < r < r max considered 
in the % 2 minimizing fit of eq. (4.1) to the lattice potential. Details on this parameter variation 
can be found in [2]. For each set of input parameters {t m \ n , t mm , r m \ n , r max ) we determine a, d 
and Eb . Then we generate histograms for a , d and Eg weighted according to the corresponding 
^ 2 /dof. The widths of these histograms are taken as systematic errors of a, d and Eb [19], while 
the statistical errors are obtained via a jackknife analysis. In Figure 1 example histograms for the 
scalar isosinglet for qq = ud are shown. 



0.28 0.3 0.32 0.34 0.36 0.38 0.4 0.42 

a 



0.2 0.3 0.4 0.5 0.6 0.7 


d in fm 


Figure 1: Flistograms for the scalar isosinglet for qq = ud. The red/green/blue bars indicate the 
statistical/systematic/combined errors. 

The resulting potentials tits for different channels, i.e. eq. (4.1) with corresponding values 
for a and d, are collected in Figure 2. The error bands represent the combined systematic and 
statistical errors. One can observe that the potentials are wider and deeper for lighter qq quark 
masses. Moreover, the scalar channels are more attractive than the respective vector channels. 
Correspondingly, it turns out that there is a bound state only for the scalar isosinglet with qq = ud 
with binding energy — Eb = 93 1 ///-. MeV, i.e. a bound state with around 2 a confidence level. 

In Figure 3 we present our results in an alternative graphical way. The three plots correspond 
to u/d, s and c light quarks qq, respectively. Each fit of eq. (4.1) to lattice potential results is 
represented by a dot (red: scalar channels; green: vector channels; crosses: r mm = 2a; boxes: r m \ n = 
3a). The extensions of the point clouds represent the systematic uncertainties with respect to a and 
d. If a point cloud is localized above or left of the isoline with Eb = —0.1 MeV (essentially the 
binding threshold), the corresponding four quarks qqbb cannot form a bound state. A localization 
below or right of that isoline is a strong indication for the existence of a tetraquark. Again there 
is clear evidence for a tetraquark state in the scalar u/d channel. The scalar s channel is slightly 
above, but rather close to the binding threshold. The scalar c and all vector channels clearly do not 
host bound four-quark states. 


5. Summary and outlook 

We have found a udbb tetraquark with quantum numbers I(J P ) = 0(1 + ) (i.e. in the scalar 
isosinglet channel with qq = ud) with a confidence level of around 2a. There seem to exist no 
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qq = uu, (ud+duyVZ, dd (vector) 


qq = ss (vector) 



r in fm r in fm 


Figure 2: Potentials fits for different channels (upper line: scalar isosinglet; lower line: vector 
isotriplet). The curves without an error band arc copied from the respective other plots in the same 
line for easy comparison. Vertical lines indicate the available lattice bb separations. 


tetraquarks for the other channels. 

In this work lattice QCD computations have been performed for light u/d quarks correspond¬ 
ing to m n ~ 340 MeV. We plan to repeat the analysis for at least another pion mass and then 
extrapolate to the physical point. It will then be most interesting to check, whether a bound state 
will also appear in the vector isotriplet channel with qq = ud. Another aspect is to investigate 
the structure of the found I(J P ) = 0(1 + ) tetraquark, i.e. to explore, whether it is rather a mesonic 
molecule or a diquark-antidiquark pair. We also plan to include corrections due to the heavy quark 
spins (for first preliminary results cf. [14]). Finally, one should study the experimentally more 
accessible, but theoretically more challenging case of qqbb systems. 
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Figure 3: Binding energy isolines Eb = constant in the a-d-plane for u/d, s and c light quarks qq 
together with the fit results of eq. (4.1) to lattice potentials. 
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